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Summary 

The results from an in-house study aimed at im- 
proving the dynamic and aerodynamic characteris- 
tics of composite rotor blades through the use of 
extension-twist elastic coupling are presented. A set 
of extension-twist-coupled composite tubular spars, 
representative of the primary load-carrying struc- 
ture within a helicopter rotor blade, was manufac- 
tured with four plies of T-650/42 graphite fiber with 
ERLX 1925-2 epoxy resin in a plain weave cloth 
prepreg. The 0°-90° cloth weave was rotated off- 
axis to achieve a [+20°] 4 laminate. These spars 
were noncircular in cross-section design and were 
therefore subject to warping deformations. Three 
different cross-sectional geometries were developed: 
D-shape, square, and flattened ellipse. Three spars 
of each type were fabricated in order to assess the 
degree of repeatability in the manufacturing process 
of extension-twist-coupled structures. Results from 
free-free vibration tests of the spars were compared 
with results from normal modes and frequency anal- 
yses of companion shell-finite-element models devel- 
oped in MSC/NASTRAN. Five global or nonshell 
modes were identified within the frequency band- 
width from 0 to 2000 Hz for each spar. Because 
of some anomalies encountered in the testing and 
subsequent determination of the modal properties 
for two of the three D-shape spars, the results for 
only one D-shape spar were determined. The fre- 
quencies and associated mode shapes were corre- 
lated with analytical results, and the agreement was 
within 13.8 percent. Frequencies corresponding to 
the five global mode shapes for the three square 
spars agreed within 9.5, 11.6, and 8.5 percent of 
their analytical counterparts. Five global modes 
were similarly identified for the three elliptical spars 
and agreed within 4.9, 7.7, and 9.6 percent of the 
analytical results. 

Introduction 

The design of advanced rotor blades is becom- 
ing significantly influenced with the incorporation 
of composite materials technology. Such designs of- 
fer considerable advantages with respect to strength 
and weight criteria, in addition to providing a means 
of efficiently controlling static and dynamic response 
through the implementation of elastic tailoring. The 
ability to analytically predict the dynamic charac- 
teristics of composite rotor blades exhibiting elastic 
coupling is essential to the implementation of elas- 
tic coupling in these applications. However, as a 
result of advancements in elastically tailored struc- 
tural design, a corresponding increase in the level of 
analytical complexity is introduced. This complex- 


ity is particularly true of the dynamic behavior, as 
the elastic couplings can have a substantial influence 
on the vibration characteristics. The identification of 
the vibration characteristics (i.e., natural frequencies 
and mode shapes) of a rotor blade is necessary for 
the prediction of its resonant behavior and is essen- 
tial in aeroelastic stability analyses. Therefore, it is 
necessary to determine the accuracy of existing anal- 
yses applicable to elastically coupled composite rotor 
blades. The validation of these analyses is further de- 
pendent upon the establishment of an extensive ex- 
perimental data base for the dynamic characteristics 
of elastically coupled structures, which has yet to be 
developed. 

Many examples of recent analytical treatments 
of composite rotor blade structures exist in the 
literature and include advanced beam and plate 
finite-element formulations. The inclusion of so- 
called nonclassical effects, which include transverse 
shear, warping, and geometric nonlinearities, has 
been emphasized in these formulations and attempts 
to account for the complex structural behavior of 
advanced helicopter rotor blades. For example, in 
reference 1, the warping effects of composite beams 
undergoing large deflections and large rotations in 
rotating and nonrotating conditions were established. 
Verification of this analysis was performed by com- 
paring numerical results with those obtained exper- 
imentally in reference 2, producing excellent corre- 
lation for flapping and lead-lag frequencies. The 
dynamic behavior of elastically coupled composite 
blades subject to large deflections was also stud- 
ied in reference 3. Several types of elastically cou- 
pled graphite-epoxy spars, including bending-twist- 
coupled and extension-twist-coupled designs for three 
varying thicknesses, were fabricated in addition to 
noncoupled designs to investigate the influence of 
elastic couplings in laminate designs. The influence 
of static deflections due to a natural gravity load on 
the small amplitude vibrations was also studied for 
each spar. Results from the experimental- analytical 
investigation showed that the presence of static de- 
flections has a significant influence on torsional and 
lead-lag modes and frequencies. This study was fur- 
ther expanded in reference 4 to investigate the effects 
of large amplitude vibrations, and similar findings 
resulted. 

The influence of ply lay-ups on the natural fre- 
quencies and mode shapes of composite beams was 
addressed in reference 5 where a mixed finite-element 
method was used. Numerical results showing the 
variation in frequency versus ply angle were pre- 
sented for a thin-walled tube with circular cross 
section for a variety of laminate designs, including 



an extension-twist-coupled laminate, where several 
variational characteristics were observed. Specifi- 
cally, some extension-twist-coupled modes occurred 
in pairs, with one mode having extension and torsion 
components of opposite signs, whereas the compo- 
nents of the companion mode shared the same sign. 
The characteristic shape of the coupled modes, as 
well as the sequential order of frequencies, was found 
to vary with ply angle. In addition, analytical re- 
sults for a laminated solid rectangular cross section 
were compared with experimental results from refer- 
ence 6 for ply lay-ups of 0°, 15°, 30°, and 90°. Good 
correlation was obtained, although only experimen- 
tal frequencies associated with the coupled vertical 
bending/torsion modes were available. 

The dynamic behavior of an extension-twist- 
coupled composite circular tube was analytically and 
experimentally studied in reference 7. An equivalent- 
beam finite-element model based on the coupled- 
beam analysis developed in reference 8 was used in 
conjunction with a detailed MSC/NASTRAN finite- 
element model to predict the natural frequencies and 
mode shapes. Experimental data from a free-free 
vibration test of a thin-walled composite test spar 
were used to validate the finite- element models. Re- 
sults from this study demonstrated that natural fre- 
quencies and mode shapes could be obtained within 
practical engineering accuracy for nonwarping-prone 
extension-twist-coupled composite structures. How- 
ever, the equivalent-beam finite-element model did 
not accurately predict the frequencies for mode 
shapes influenced by elastic coupling. 

The global response characteristics of general 
thin-walled anisotropic beams which exhibit shear 
and warping deformations were studied in refer- 
ence 9. Results from an anisotropic beam finite- 
element analysis were experimentally verified with 
a set of static axial tests performed on differ- 
ent configurations of extension-twist-coupled tubular 
beam structures to assess the influence of warping- 
prone cross-sectional geometries. Test results demon- 
strated that bending-related warping effects have sig- 
nificant influence on beam-like structures and are 
amplified by the extension-twist coupling, whereas 
the effects of torsion-related warping were not found 
to be significant for practical rotor structures. This 
investigation, however, only addressed the static be- 
havior of extension-twist-coupled composite spar de- 
signs and leaves the question of how such designs 
influence the dynamic characteristics. 

An in-house research program aimed at im- 
proving the dynamic and aerodynamic character- 
istics of composite rotor blades through the use 
of elastic couplings has been initiated to address 


this need. Extension- twist elastic coupling, which 
exhibits coupling among the extensional and tor- 
sional stiffnesses of a structure, was chosen for this 
study because of its potentially advantageous appli- 
cation to tilt rotor blades. The objectives of this 
program were to provide a basis for understanding 
the mechanics of elastic coupling as applied to rotor 
blade structures and also to develop verifiable finite- 
element modeling methodologies that are capable of 
accurately representing this phenomenon. A series 
of analytical and experimental studies has been de- 
signed to demonstrate the potential improvements 
to tilt rotor blade performance. Analytical models 
would be validated through the static and dynamic 
testing of a series of elastically coupled composite test 
spars representative of rotor blade spars. These test 
spars include untwisted, thin- walled tubes of circular 
cross section and progress to more advanced designs 
involving noncircular cross sections and pretwist. 
Such designs introduce additional complexity as they 
are subject to out-of-plane warping deformations. 
This will lead to the development of low twist, model- 
scale composite rotor blades, with subsequent whirl 
and wind-tunnel tests designed to demonstrate the 
feasibility of changing blade twist as a function of 
rotor speed. The analytical-experimental study of 
these blades will provide the basis for the design and 
development of a highly twisted composite tilt ro- 
tor blade design incorporating extension-twist elastic 
coupling. 

This paper presents recent results from an in- 
house investigation of the dynamic characteristics of 
thin-walled extension-twist-coupled composite tubu- 
lar spars. Results from a normal modes analysis of 
shell-finite-element models were compared with ex- 
perimental results obtained from free-free vibration 
tests of extension-twist-coupled composite tubular 
spars. These spars were noncircular . in cross-section 
design and were therefore subject to warping defor- 
mations. Because the presence of extension-twist 
coupling in a laminate acts to “magnify” the warp- 
ing effect, a significant influence on the dynamic be- 
havior may result. Such warping-prone designs are 
more closely representative of true rotor blade spar 
structures and require thorough investigation of the 
dynamic characteristics. 

Symbols 

Eh extension modulus (spanwise), lb/in 2 

E 22 extension modulus (chordwise), lb/in 2 

FEA finite-element analysis 

G \2 shear modulus, lb/in 2 
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j 

response degree of freedom 

MAC 

modal assurance criteria 

MCF 

modal confidence factor 

MMIF 

multivariate mode indicator function 

m 

number of response degrees of freedom 

*ply 

lamina thickness, in. 

1/12 

Poisson ratio 

P 

density, lb/in 3 


771 -dimensional vector representing 
normalized mode shape 1 

V>2 

ra-dimensional vector representing 
normalized mode shape 2 

Superscript: 


* 

complex conjugate 

Test Models 


Composite Test Spars 

Nine composite tubular spars were fabricated for 
dynamic testing. Each tube was 24 in. in length and 
was fabricated with four plies of T-650/42 graphite 
fiber with ERLX 1925-2 epoxy resin in a plain weave 
cloth prepreg. This material system is manufactured 
by Amoco Performance Products, Inc. The associ- 
ated material properties, obtained from coupon tests 
conducted in reference 9, are shown in table 1. The 
0°-90° cloth weave was rotated off-axis to achieve 
a [+20°]4 ply stacking sequence. This symmetric 
arrangement of ply angles provides a mechanism 
for extension-twist coupling in a thin-walled cellu- 
lar structure through the interaction of extensional 
and torsional deformations. Specifically, the lami- 
nate shears in response to an axial load and causes 
the cell structure to develop a net twist deformation. 
In addition, a secondary type of coupling that was in- 
herent in this type of extension-twist-coupled design 
acted to couple the bending and shear deformations. 
This parasitic bending-shear coupling influences a 
structure such that a shear deformation in a given 
geometry plane results from a bending deformation 
in an orthogonal plane. 

Three different cross-sectional geometries were 
manufactured: D-shape, square, and flattened el- 
lipse. In order to assess the degree of repeatability in 
the manufacturing process of extension-twist-coupled 
structures, three spars of each type cross section were 
fabricated. These spars were given a three letter des- 
ignation Vci, where V denotes vibration test arti- 
cle, c refers to cross section (i.e., D-shape, Square, 


Elliptical), and i refers to spar number (i.e., 1, 2, 3). 
Figure 1 depicts the three cross-sectional shapes with 
the associated dimensions. Although all spars were 
designed to have a nominal wall thickness of 0.030 in., 
the actual wall thicknesses varied slightly because of 
the manufacturing process. Because an exact span- 
wise variation of the wall thickness was difficult to 
accurately measure by nondestructive means, an av- 
erage wall thickness was computed for each spar by 
dividing the weight of the spar by the product of 
the associated material density, spanwise length, and 
contour length. This average thickness was used in 
the finite-element analysis models of the spars. A 
summary of the average wall thicknesses determined 
in this manner as well as the percent difference from 
the design value appears in table 2 along with the 
measured weight for each spar. The deviation in wall 
thickness was generally quite small (within 2.3 per- 
cent of the design value), except for the elliptical 
spars. Average deviations of approximately 10 per- 
cent were observed in this subset of spars. Poten- 
tial causes for the marked differences in the elliptical 
spar wall thicknesses may have precipitated from a 
change in fabrication personnel during this part of 
the manufacturing phase. Factors such as the tem- 
perature rise rate during the cure cycle and wrap- 
ping of the prepreg during the lay-up procedure are 
integral in the fabrication process and inconsisten- 
cies in either or both of these variables could have 
conceivably contributed to the differences noted. 

Spar Fabrication Process 

The thin-walled tubular spars were fabricated by 
utilizing an expandable rubber mandrel and a seg- 
mented aluminum female mold configured to the de- 
sired external spar geometry. The rubber mandrel 
was sized to accommodate the composite wall thick- 
ness and a free lateral expansion gap of approxi- 
mately 0.020 in. prior to contact with the surface 
of the female mold during processing. A (mild) steel 
rod was placed in the center of the mandrel to provide 
stiffness for handling purposes after consolidation of 
the room temperature vulcanizing silicon rubber. Af- 
ter the rubber was cast, cured, and removed from the 
sizing tool, it was sprayed with a wet coat of fluoro- 
carbon release agent and dried in an air circulating 
oven. The mandrel was then ready for application of 
the graphite- epoxy cloth composite. 

Figure 2 shows the silicon rubber mandrel and the 
initial wrap of composite material which has been 
cut to the desired fiber orientation. The wrapping 
procedure began by aligning the rubber tool with 
one edge of the precut composite material as shown 
in figure 3. Subsequent layers of composite material 
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were applied to the assembly by continuing to wrap 
successive layers of composite in the proper fiber 
orientation, as shown in figure 4. Once the desired 
thickness was built up onto the rubber mandrel, the 
assembly was placed in the female aluminum tool 
cavity (fig. 5) and the cover plate was secured by 
bolts. A control thermocouple was inserted into the 
material at one end and the tool was placed in an 
air circulating oven. The temperature of the oven 
was increased at a rate of approximately 10°F/min 
until the spar reached the cure temperature of 350°F. 
Upon heating, the rubber mandrel expanded, which 
forced the composite material against the surface of 
the female mold. The rubber maintained a pressure 
of approximately 200 psi on the composite laminate 
at 350°F. This condition was held for 1 hour for 
complete composite cure. The assembly was cooled, 
removed from the oven, and disassembled. The 
rubber core shrank during the cooling process, which 
allowed free removal from the composite tube spar. 
Examples of the cured composite tube spars are 
shown in figure 6. 

Richard Bennett, Materials Processing and Fab- 
rication Division, participated in the manufacturing 
of the elliptical spars. 

Test Results 

Vibration Test Setup 

The vibration test setup used to determine the 
dynamic characteristics of the composite tubes is 
shown in figure 7. The test spars were suspended 
from a support structure with light rubber tubing 
to simulate a free- free boundary condition. A cal- 
ibrated impact hammer provided the input excita- 
tion for the spar. The test procedure consisted 
of a roving hammer concept such that impacts at 
multiple locations along the structure were made 
while acquiring response data from a fixed loca- 
tion. A three-axis accelerometer transducer config- 
uration was mounted on one end of each spar. (See 
fig. 8.) The configuration incorporated three single- 
axis piezoelectric accelerometers, which were affixed 
to an aluminum mounting block to form a triaxial 
accelerometer arrangement. This served as the ref- 
erence point for subsequent measurements. Smaller 
blocks were used on the elliptical and D-shape spars 
to facilitate vertical and lateral impacts on the flat 
and curved surfaces. A list of the block types and 
their corresponding weights is given in table 3. 

A prescribed set of excitation points was chosen to 
coincide, as closely as possible, with the correspond- 
ing grid points of the finite-element model which was 
formed for each spar. By doing this, an accurate 


representation of grid point deformations could be 
obtained when comparing test and analysis models. 
The excitation and response point degrees of freedom 
(location and direction) for each spar are listed nu- 
merically in tables 4 through 10 and are illustrated 
for each spar type in figure 9. Twelve locations were 
excited in both the lateral (F) and vertical (Z) di- 
rections for each square spar by using a configura- 
tion of two parallel rows of six axial locations. A 
similar arrangement was employed for the D-shape 
and elliptical spars, although vertical ( Z ) inputs were 
omitted along the leading edge of the D-shape spar 
(grid points 7-12), as well as the lateral (F) inputs 
along the top of the elliptical spars (grid points 1-6). 
In the latter case, significant difficulty was encoun- 
tered when attempting lateral impacts on the up- 
per (flat) surface of the ellipse because of the low 
shearing modulus of the cyanoacrylate adhesive used 
to mount the aluminum blocks. Therefore, only im- 
pacts which were perpendicular to the shearing plane 
of the adhesive were obtained at these locations. Ad- 
ditionally, a slight shift in excitation point location 
from the desired location was necessary in a few in- 
stances to achieve a cleaner frequency input spec- 
trum for a spar. This caused a slight nonuniformity 
in excitation point locations. This shift was usually 
required whenever slight anomalies in the composite 
skin were encountered (i.e., excessive accumulation 
or depletion of resin, inconsistent fabric weave). 

Data Acquisition System 

The vibration test of the spars was performed 
with the Hewlett-Packard 3565S modal testing/data 
acquisition system. (See fig. 10.) This system is 
self-contained with respect to hardware components 
and consists of a front end (i.e., input and source 
modules, charge amplifiers, signal processing mod- 
ules), central processing unit (cpu), fixed disk drives, 
an optical disk drive, and a magnetic tape drive. 
The front end is capable of simultaneously acquir- 
ing up to 8 channels of data with a frequency range 
of 50 kHz. 

Modal data acquisition was performed by the 
front end, under the control of the Hewlett-Packard 
Virtual Instrument for Test and Analysis (HP- 
VISTA) signal processing software (ref. 10). The 
VISTA software was used to collect data from the 
measurement hardware and served as a tool for anal- 
ysis of measurement results. It effectively trans- 
formed the workstation into an “analyzer,” or virtual 
instrument. The raw analog signals generated by the 
accelerometers or force transducers were fed into the 
front end and passed through analog-to-digital filters 
and antialiasing filters. The resulting digital signals 
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were then processed by the fast Fourier transform 
(FFT) algorithm and transformed into the frequency 
domain. The transfer functions were stored as binary 
files on disk for subsequent processing and analysis 
and/or displayed on the monitor for viewing. 

Test data were obtained for each spar over the 
frequency bandwidth from 0 to 3200 Hz. The block 
size, which determined the number of digitized sam- 
ples within the time record, was 4096 (real), with a 
frequency resolution of 2.0 Hz. The associated time 
length of actual data acquisition was 500 msec/data 
average. A total of 30 data averages was obtained at 
each excitation point location for all spars. 

The associated data files were then ported to the 
modal analysis software, I-DEAS Tdas (test data 
analysis software), which is a general purpose test 
data analysis package developed by Structural Dy- 
namics Research Corporation (ref. 11). The porting 
of these data was facilitated by a Hewlett-Packard 
Modal Data Manager, a utility software program 
(ref. 12), which assembled the individual frequency 
response function files created by VISTA to produce 
a formatted file in which a matrix of frequency re- 
sponse functions was stored. This file format was 
the means by which data were transferred to I-DEAS 
Tdas, which extracted modal parameters (frequency, 
mode shape, and damping information) from the 
transfer functions. 

Data Reduction Techniques 

The data reduction phase of this study is per- 
haps most easily explained in terms of five proce- 
dure steps involving curve-fitting techniques, esti- 
mating curve-fit parameters, determining residues, 
calculating mode shapes, and correlation. 

Step 1 — curve- fitting techniques . The Tdas 
software provided several curve-fitting techniques for 
the extraction of frequency, amplitude, damping, and 
phase information for a given mode. These tech- 
niques included single-degree-of-freedom curve fits, 
such as circle fit and peak or valley searches, as well 
as multiple-degree-of-freedom curve fits, which in- 
cluded complex exponential, direct parameter, and 
poly reference techniques. The poly reference curve- 
fitting technique (ref. 13) was chosen to extract the 
modal parameters for the test spars. This time do- 
main curve-fitting technique uses the impulse re- 
sponse function (inverse of the Fourier transform 
of the frequency response function) from single or 
multiple excitation/reference locations in a global 
least-squares manner. This technique is advanta- 
geous because all the available frequency response 
functions may be used to calculate a global esti- 
mate for frequency, damping, and residue; thereby 


a more complete description of the modal system is 
provided. The damped natural frequencies and cor- 
responding damping ratios are first extracted from 
the acquired time data, followed by calculation of 
the residues (amplitude and phase). The polyrefer- 
ence technique is useful for multiple-input/multiple- 
output (MIMO) applications and situations where 
closely coupled or coincident modes exist. The curve- 
fitting procedure for the polyreference technique is 
an iterative process, involving the construction of a 
correlation matrix which identifies the reference co- 
ordinates to be used (i.e., grid point number and di- 
rection), the analysis bandwidth, and the maximum 
number of roots allowed in the estimation. The ac- 
curacy of the polyreference technique is dependent 
upon user-specified estimates of the number of roots 
present in the frequency range of interest. 

Step 2 — estimating curve-fit parameters . 
An analytical tool designed to aid the user in the de- 
termination of the root estimate, known as the multi- 
variate mode indicator function (MMIF), was used 
in this study. The MMIF can detect the presence of 
modes within the frequency range of interest, as well 
as repeated roots, by using the frequency response 
functions from a given set of reference (response) lo- 
cations. Conceptually, the MMIF can be thought of 
as a graphical representation of the ratio of the ki- 
netic energy of in-phase responses to the kinetic en- 
ergy of the total response plotted as a function of fre- 
quency. This ratio is minimized so that the real part 
of the response is forced to be as small as possible 
compared with the total response. Minima observed 
in this plot indicate real normal modes that can be 
excited within the structure for a given set of refer- 
ence locations. The number of mode indicator func- 
tions in the MMIF is equal to the number of reference 
coordinates chosen, up to a maximum of three. Three 
reference coordinates, 1X+, 1Y_, \Z+ (grid point 1 
in the positive X , negative Y y and positive Z di- 
rections, respectively), were selected in this study; 
therefore, three functions (primary, secondary, and 
tertiary) were calculated. The MMIF can detect re- 
peated roots which may exist in the structure via co- 
incident minima in the secondary and tertiary func- 
tions. A more detailed derivation and explanation of 
the MMIF can be found in reference 14. 

Once the MMIF has been calculated, the num- 
ber of minima as well as the corresponding natu- 
ral frequencies can be determined. This determina- 
tion, which can be performed by visual inspection or 
through a curve-fit procedure, establishes an initial 
estimate for the number of roots within the frequency 
range for the polyreference curve-fit method. The 
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estimate of roots is an iterative process and judg- 
ment is required by the analyst because computa- 
tional roots (roots which have no physical signifi- 
cance) are also produced by the curve fit and must 
be identified and eliminated. 

Step 3 — determining residues . Once the fre- 
quency and damping terms have been estimated, real 
or complex residues (amplitude and phase) are calcu- 
lated based on the specified references and responses. 
Residues of the real type imply that the phase angles 
of the associated mode shape differ by 0° or 180° such 
that all the structural points physically reach their 
maximum displacements at the same time. Real- 
type modes are based on proportional damping. A 
mode shape of complex type exhibits nonpropor- 
tional damping, where the phase angles can have 
values other than 0° or 180°, and therefore physi- 
cally reach their maximum displacements at different 
times. In instances where damping is light, such as 
that anticipated in this investigation, the assumption 
of real- type modes is generally valid. 

A feature within the polyreference algorithm 
known as the modal confidence factor or MCF 
(ref. 15) was used as an aid to discriminate between 
computational and physical modes. The MCF takes 
advantage of the redundant phase relationships that 
are satisfied by physical modes but not by compu- 
tational modes and, therefore, represents a measure 
of how well a given mode conforms to this relation- 
ship. The MCF assumes a value near unity for phys- 
ical modes and a value less than unity (typically 
less than 0.8) for computational modes. The theory 
behind the MCF is derived in reference 15. 

Step 4 — calculating mode shapes . Mode 
shapes corresponding to the real roots identified in 
the previous step were generated and scaled to a 
user-specified response coordinate, usually a driving 
point measurement. The mode shapes for all spars 
were scaled to their respective 1 driving point re- 
sponse, which is defined as the measurement made 
at the same degree of freedom as that of the excita- 
tion (i.e., grid point 1 in the negative Z direction). 
The resultant mode shapes were displayed and ani- 
mated with simplified graphics models representative 
of each of the test spars. These models consisted of 
grid points which geographically corresponded to the 
excitation points on the spar and were connected by 
plate elements as a means of defining the spar out- 
line. Because data were obtained at only two lateral 
locations for each spar (fig. 9), the graphics models 
appeared planar (two-dimensional). 

Step 5 — correlation . Once the modal proper- 
ties of the structures were determined, a comparison 


of the test and analysis mode shapes was performed 
with the modal assurance criteria (MAC) matrix. 
The MAC matrix, which is a three-dimensional ma- 
trix, is useful for identifying mode shapes that are 
similar to each other and can be calculated for test 
versus test, analysis versus analysis, or test versus 
analysis mode shapes. The matrix is constructed 
such that the rows and columns of the matrix cor- 
respond to test and/or analysis mode shapes, and 
the vertical axis corresponds to the associated MAC 
value. The MAC value is a scalar ranging from zero 
to unity and represents the linear dependence or sim- 
ilarity of a pair of mode shapes, with unity indicating 
perfect similarity and zero indicating total dissimilar- 
ity. Although this value may be regarded as a mea- 
sure of orthogonality, it does not constitute a true or- 
thogonality check, because the associated mass and 
stiffness matrices are not considered in the calcula- 
tion of the MAC value. The MAC value is calculated 
by the following equation from reference 15: 


MAC = 
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Although the MAC matrix is a very useful tool in 
performing mode shape correlation, user judgment 
is required, as contaminating effects such as geomet- 
ric aliasing, nonlinearities, or noise may contribute 
to the calculation of erroneous MAC values. As 
pointed out in reference 16, the MAC value is only 
a measure of consistency, and not necessarily valida- 
tion, between modal vector estimates. This is be- 
cause the MAC cannot account for errors that may 
be consistent among all modal vector estimates. 

Analysis Models 

The composite spars were modeled and analyzed 
with the MSC/NASTRAN finite-element program 
(ref. 17). Because each spar varied slightly with re- 
spect to the nominal wall thickness and weight pre- 
scribed for the spar type, a separate finite-element 
model was formed for each spar. With this ap- 
proach, it was possible to more closely model each 
spar and account for individual variations than with 
a single model which has the averaged or “smeared” 
properties for each spar type. The finite-element 
models for each spar type are shown in figures 11 
through 13. The structural modeling of each spar 
consisted of two-dimensional, quadrilateral flat plate 
(CQUAD4) elements having anisotropic properties. 
Although separate models were constructed for each 
spar, the same number of elements and grid points 
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was used for each model of a given spar type. All 
models were formed with 25 plate elements in the 
lengthwise direction, with the number of circumfer- 
ential elements varying with cross-section type. (See 
appendix for results of a limited convergence study.) 
Eight plate elements were used circumferentially to 
model the square cross section, whereas 12 and 14 el- 
ements were used for the D-shape and elliptical ge- 
ometries, respectively. The total degrees of freedom 
for each model were reduced through static conden- 
sation, with a reduced set of dynamic degrees of free- 
dom resulting. The total number of grid points, as 
well as the number of plate elements and dynamic 
degrees of freedom for each model, is also listed in 
the associated figures. The plane of grid points for 
each model was located at the geometric midplane 
of the 4-ply laminate and was therefore coincident 
with the element reference plane. Structural mass 
was accounted for through the use of the consistent 
mass matrix option in the analyses. The nonstruc- 
tural mass and associated inertia of the accelerome- 
ter mounting blocks attached to each spar were ac- 
counted for through the use of concentrated mass 
elements and solid elements. Concentrated mass 
(C0NM2) elements were used to model the small 
impact blocks used for excitation input on the ellip- 
tical and D-shape spars, which used two blocks and 
one block, respectively. The accelerometer mount- 
ing block used on the square and D-shape spars 
was modeled with a single (CHEXA) solid element. 
The mounting block used on the elliptical spars was 
modeled by using four (CPENTA) solid elements. 

Results 

Experimental Results 

D-shape spar . Although three D-shape spars 
(VD1, VD2, and VD3) were fabricated for this in- 
vestigation, some anomalies were encountered in the 
testing and modal analyses for two of the three spars 
(VD2 and VD3). The cause of these anomalies, 
which were manifested as multiply-split poles and fre- 
quency shifts, could not be explained satisfactorily, 
although it may be likely due to material or lamina- 
tion defects within the spars in question. Therefore, 
the two suspect spars were not considered further in 
the investigation and results are not presented for 
these spars. 

The multivariate mode indicator function for the 
VD1 spar was calculated over the frequency range 
from 0 to 2000 Hz, where the driving point coherence 
exhibited the best characteristics and is shown in fig- 
ure 14. Considering the relative magnitudes of the 
primary minima, approximately 15 modes appeared 


to be within the range of interest. Although there is 
no threshold value which explicitly defines a minima 
as a potential mode, those minima occurring within 
the lower third of the value scale may generally be 
considered valid in this example. Minima occurring 
near this limit and above may also be considered al- 
though user judgment is required within this “gray 
area” in light of the rather ragged quality of the pri- 
mary function, particularly in the frequency range 
from 0 to 300 Hz. This raggedness is most likely due 
to low coherence values which were observed in the 
cross measurements involving the axial (x) measure- 
ments in this range. Apparent double- peaked min- 
ima occurring in the mid to upper ranges of the scale, 
as well as the slight drops in the secondary function, 
may likely be influenced by slight differences or shifts 
in frequency which were sometimes observed between 
measurement locations. No significant evidence of re- 
peated roots was found, which would be manifested 
by substantial and coincident drops in the secondary 
and possibly tertiary functions in the graph. 

With the information obtained from the MMIF, 
the analytical curve fit of the experimental data 
was performed by using the polyreference technique. 
The resulting synthesized frequency response func- 
tion, plotted as an overlay on the experimental fre- 
quency response function, is shown in figure 15. The 
analytically synthesized function is nearly identical 
to its experimental counterpart throughout the fre- 
quency range from 0 to 2000 Hz, with the excep- 
tion of the mode or modes near 1950 Hz. The phase 
has deteriorated considerably in this frequency range; 
this results in diminished accuracy in the analytical 
function. The corresponding frequency and mode 
shape data extracted from the polyreference curve 
fit in the real domain are listed in table 11, including 
the estimates of damping, phase angle, amplitude, 
and MCF. (Modal parameters were also obtained 
with the polyreference curve fit in the complex do- 
main and did not differ significantly from those ob- 
tained in the real domain.) Thirteen modes were ex- 
tracted within the frequency range from 0 to 2000 Hz. 
(Although the MMIF indicated potential modes oc- 
curring near 1850, 1925, and 1980 Hz, these min- 
ima were determined to correspond to local and/or 
frequency-shifted data and could not be captured by 
curve fits in a global or reduced bandwidth.) Of 
these 13 modes, 5 global “nonshell” modes, which 
are the primary mode types of interest, were identi- 
fied: first vertical bending, first lateral bending, sec- 
ond vertical bending, first torsion, and second lateral 
bending. The mode shapes were scaled to the 1 
driving point measurement. Because a frequency 
response function is a complex quantity (frequency 
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plus a real and an imaginary part), a set of mul- 
tiple plots is required to completely represent the 
frequency response function because to fully de- 
scribe it is not possible when using a single two- 
dimensional plot. One method of presenting the fre- 
quency response data involves plotting the real and 
the imaginary parts of the function versus frequency. 
This type of plot illustrates how the phase changes 
through a resonance, which is defined as the point 
where the real part is zero and the imaginary part has 
a maximum value. The magnitude of the frequency 
response function is equivalent to the square root of 
the sum of the squares of the real and imaginary 
parts. An additional plot associated with frequency 
response measurements is the coherence, which can 
be thought of as a ratio of signal to noise. The 
coherence represents a measure of the causality be- 
tween the input and response signals and is always 
less than or equal to unity. A more in-depth discus- 
sion of these plots, including the mathematical basis, 
can be found in reference 18. Plots of the 1 Z- driv- 
ing point measurement are shown in figure 16 using 
the plot types described above. The global mode 
shapes are plotted in figure 17, with the dashed out- 
line representing the undeformed shape. Because the 
grid points along the leading edge of the D-shape 
spar were not impacted in the vertical (Z) direction, 
they appear undeformed in the plotted mode shapes. 
This is particularly evident in the mode shapes which 
have primarily vertical bending motions. Therefore, 
this incongruity must be considered when correlat- 
ing the test and analysis mode shapes. The bending 
deformations of the D-shape spar appear slightly cou- 
pled with respect to the vertical ( X-Z ) and lateral 
( X-Y ) planes. This coupling is a result of the para- 
sitic bending-shear coupling which is inherent in an 
extension-twist-coupled laminate design. Bending- 
shear coupling influences a structure such that a 
shear deformation in a given plane is due to a bend- 
ing deformation in an orthogonal plane. With the 
first lateral bending case (fig. 17(b)) as an example, 
a slight (secondary) deformation in the X-Z plane 
is observed in response to the primary bending de- 
formation in the X-Y plane. This effect is observed 
in each of the bending mode cases, where the sec- 
ondary deformation appears to assume a shape char- 
acteristic of the next higher bending mode (i.e., a 
slight second bending-type deformation is observed 
in the secondary plane in response to a first bending- 
type deformation in the primary plane). Also, the 
node points (points of zero displacement) are slightly 
asymmetric with respect to the center of the beam. 
This is because the application of the aluminum block 
at one end of the spar has shifted the beam mass 
center. Since this promotes a global change in the 


beam mass distribution, this asymmetry will there- 
fore exist in all observed modes. The shape of the 
first torsion mode (fig. 17(d)) is masked by the lack 
of participation of the leading-edge grid points in the 
vertical direction; therefore the determination of the 
mode shape is complicated. A torsion deformation is 
apparent in the vertical (X-Z) plane with the change 
of phase between the right and left halves of the spar 
(with respect to the center). The lateral (X-Y) plane 
perspective shows a slight lateral bending deforma- 
tion over the grid points closest to the coordinate 
system origin. In this view, the mode shape assumes 
more of a vertical shearing characteristic than a tor- 
sional one. 

Square spars . The experimental results for the 
square tubes (VSl, VS2, and VS3) are presented in 
this section. The multivariate mode indicator func- 
tion is presented for each of the three square spars 
in figures 18, 19, and 20, respectively. In each case, 
the MMIF plots exhibited very similar characteris- 
tics and indicated approximately 10 modes within 
the frequency range from 0 to 2000 Hz for VSl and 
VS3, although an 11th shell-type mode was apparent 
for VS2. Although several closely spaced modes were 
apparent, there was again no evidence of repeated or 
coincident roots. The synthesized frequency response 
function for VSl is plotted with the associated ex- 
perimentally derived function in figure 21. Overall, 
a very good fit was obtained throughout the range 
from 0 to 2000 Hz, with peaks as well as valleys ac- 
curately reproduced. Ten modes were extracted for 
VSl and are listed in table 12. The mode shapes 
were scaled to the 1 Z- driving point measurement, 
which is shown in figure 22. Five global modes were 
identified for VSl: first vertical bending, first lat- 
eral bending, second vertical bending, second lateral 
bending, and first torsion. The experimental mode 
shapes are illustrated in figure 23. The first vertical 
and lateral bending modes appear slightly coupled 
with their respective orthogonal planes such that the 
first mode bending deformation was apparent in both 
planes because the cross-sectional axes about which 
bending occurred were canted roughly 20° from the 
true vertical and lateral axes, respectively. The iner- 
tial asymmetry produced by the accelerometer block 
is fostering this effect, as the bending deformations 
would be otherwise uncoupled. (This was analyti- 
cally confirmed by removing the block in the finite- 
element model and recalculating the mode shapes. 
No evidence of this coupling action was apparent in 
the recalculated mode shapes.) The first vertical and 
lateral bending modes appear uncoupled with respect 
to torsional deformation. However, it is noted that 
some torsional deformation is present in each of the 
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second bending modes. This coupled behavior is be- 
lieved to be a function of the aforementioned iner- 
tial asymmetry because the spar is not structurally 
bending-twist coupled by design. The first torsion 
mode was similarly influenced and exhibited a slight 
second bending-type deformation in the lateral plane. 

The analytically synthesized frequency response 
function is plotted as an overlay on the experimental 
frequency response function for the VS2 spar in 
figure 24. Eleven modes were extracted for the 
VS2 spar, as listed in table 13. The mode shapes 
were scaled to the 1 Z- driving point measurement, 
which is shown in figure 25. Five global modes were 
identified: four bending modes and one torsion mode. 
The experimental mode shapes are nearly the same 
as those of VSl and are thus not shown. 

The experimental and analytically synthesized 
frequency response functions are shown in figure 26 
for the final square spar, VS3. Ten modes were 
extracted within the frequency test range from 0 
to 2000 Hz and are listed in table 14. These mode 
shapes were also scaled to the 1 Z- driving point mea- 
surement shown in figure 27. As with the other two 
square spars, four global bending modes and the tor- 
sion mode were identified for VS3. The experimen- 
tal mode shapes for VS3 were again similar to those 
obtained from VSl and are not shown. 

Elliptical spars. The test results for the three 
elliptical spars (VE1, VE2, and VE3) are presented 
in this section. The corresponding multivariate 
mode indicator functions for VE1, VE2, and VE3 
are shown in figures 28, 29, and 30, respectively. 
Seven modes were evident within the frequency range 
from 0 to 2000 Hz for VE1, whereas six modes were 
evident for VE2 and VE3. Note the coincident drops 
in the secondary MMIF functions for VE2 and VE3, 
which appear to indicate coincident or repeated roots 
at the associated frequencies, given the somewhat 
significant magnitude of the minima. 

The corresponding analytically synthesized fre- 
quency response function, as obtained from the 
polyreference curve fit, is plotted with the experi- 
mental frequency response function in figure 31 for 
the VE1 spar. In general, the curve fit exhibited 
very good characteristics within the frequency range 
from 0 to 2000 Hz although the antiresonance near 
1750 Hz appeared slightly shifted between the test 
data and the analytical curve fit. Seven modes were 
determined in the frequency range from 0 to 2000 Hz 
for the elliptical spar and are listed in table 15. Al- 
though the MMIF function, shown in figure 28, indi- 
cated a potential occurrence of repeated roots within 
the frequency range, none were found. The mode 


shapes were scaled to the 1 Z- driving point response 
shown in figure 32. A significant deterioration in the 
coherence function was evident near 1600 Hz for the 
driving point response, perhaps affecting the qual- 
ity of the curve fit in this region. Five global modes 
were identified for VE1 within the frequency range 
from 0 to 2000 Hz: first vertical bending, first lat- 
eral bending, second vertical bending, second lateral 
bending, and first torsion. The experimental mode 
shapes for VE1 are shown in figure 33. The bend- 
ing modes are classical in nature and are not coupled 
to torsional deformation. The effects of the para- 
sitic bending-shear coupling are again evident in the 
bending modes. The grid points lying on the up- 
per surface of the elliptical spar (grid points 1-6) 
did not receive a lateral excitation and therefore ap- 
pear undeformed. The first torsion mode exhibits 
a slight coupling motion in the vertical ( Z ) plane 
such that a near second bending deformation is man- 
ifested. A sixth global mode characteristic of third 
vertical bending was strongly evident in the mode at 
1865.52 Hz, shown in figure 33(f). 

The synthesized and experimental frequency re- 
sponse functions for the VE2 spar are plotted in 
figure 34. The overall agreement is very good in 
the frequency range from 0 to 2000 Hz. Six modes 
were extracted from the poly reference curve fit, with 
no evidence of repeated roots, and are listed in ta- 
ble 16. Five global modes were identified, with simi- 
lar shape and modal characteristics as compared with 
VE1, and were scaled to the 1 Z_ driving point mea- 
surement shown in figure 35. A sixth global mode 
characteristic of third vertical bending was evident 
in the mode corresponding to 1739.82 Hz. The ex- 
perimental mode shapes are not shown because of 
their similarity with those of VE1. 

The analytically synthesized and experimental 
frequency response functions are shown in figure 36 
for the VE3 spar. Once again, the magnitude and 
phase is very good throughout the frequency range 
from 0 to 2000 Hz, with excellent representation of 
all peaks and phase changes. The modal proper- 
ties were subsequently determined for VE3, where 
six modes were determined in the frequency range 
from 0 to 2000 Hz. The resulting parameters for 
the six modes are listed in table 17. The 1 Z_ driv- 
ing point measurement for VE3 was used to scale 
the mode shapes and is shown in figure 37. Five 
global modes, similar to those obtained for VE1 and 
VE2, were identified. Additionally, the mode at 
1838.58 Hz appeared to be characteristic of third 
vertical bending. 
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Analytical Results 

D-shape spar . The results from the nor- 

mal modes and natural frequency analysis of the 
MSC/NASTRAN finite-element model of VD1 are 
listed in table 18 for the frequency range from 0 
to 2000 Hz. The calculated frequencies which cor- 
respond to the global modes identified in the exper- 
iment are indicated in table 18. The corresponding 
mode shapes are illustrated in figure 38. Deforma- 
tion of the cross sections adjacent to the accelerome- 
ter block is apparent for most of the analytical mode 
shapes. The magnitude of this deformation is greatly 
exaggerated in the figure because of the scaling of the 
plot. Also, the solid element appears to be stretch- 
ing or deforming in the global mode shapes, but 
this is again an artifact from the plot scaling. Lo- 
cal boundary effects between the block and the plate 
element may further contribute to this distorted de- 
formation. In order to more accurately capture the 
specific cross-sectional deformation, an increase in 
the number of plate elements around the circumfer- 
ence of the spar would be necessary. The calculated 
global bending modes are similar in shape to those 
obtained experimentally. The second vertical and 
lateral bending modes exhibit a slight torsional defor- 
mation, particularly evident along the leading edge 
of the D-shape, which is not apparent in the first 
bending modes. The torsional characteristic in the 
mode at 1254.97 Hz is very evident in figure 38(d), 
where most of the deformation occurs along the lead- 
ing edge. This is expected based on an anticipated 
increase in the structural stiffness near the corners of 
the D-shape. Conversely, the torsional deformation is 
slightly coupled to a second vertical bending-type de- 
formation, especially when viewed in the X-Z plane. 

Square spars . The frequencies obtained from 
a normal modes analysis for the square spars (VSl, 
VS2, and VS3) are shown in table 19 for the fre- 
quency range from 0 to 2000 Hz. The frequencies 
corresponding to global modes are indicated in the 
table. The mode shapes for VSl are depicted in 
figure 39. 

The first vertical and lateral bending modes of the 
square spars are classical in nature and do not appear 
to be significantly affected by the inertia of the ac- 
celerometer block. The second vertical ( Z ) bending 
mode is only slightly affected by the associated cross- 
sectional deformation near the accelerometer block 
and this deformation is characteristic of a vertical 
shearing. In addition, the second bending modes 
are coupled to torsional deformation, perhaps as a 
result of the inertial asymmetry produced by the ac- 
celerometer block. The first torsion mode is apparent 


at 1753.94 Hz and appears slightly coupled to second 
bending deformation. 

The resultant analytical mode shapes for the VS2 
and VS3 spars are nearly the same as those obtained 
for VSl and are thus not shown. 

Elliptical spars . The frequency and mode shape 
results for the elliptical spars (VE1, VE2, and VE3) 
are listed in table 20. The global mode shapes are 
indicated. The calculated mode shapes for VE1 are 
shown in figure 40. 

The first vertical and lateral bending modes ex- 
hibit classical bending behavior, with little cross- 
sectional deformation. The same is generally true 
of the second vertical and lateral bending modes,' al- 
though these modes appear to be slightly coupled 
with torsional deformation; this could be due to the 
asymmetrically mounted accelerometer blocks. In 
addition, the first two cross sections in the analyt- 
ical model appear somewhat flattened, although the 
magnitude of this deformation is a function of scaling 
and is indeterminate. The first torsion mode, shown 
in figure 40(e), exhibits a slight coupling with a sec- 
ond bending- type deformation. The exaggerated size 
of the cross section is again a function of the scaling. 

The corresponding mode shapes calculated for the 
VE2 and VE3 spars are nearly the same as those 
obtained for VE1 and are thus not shown. 

Experimental and Analytical Correlation 

As a means of facilitating test-analysis correla- 
tion, the analytical model geometries generated in 
MSC/NASTRAN for each spar were cross mapped 
onto the corresponding experimental model geome- 
tries. Grid points common to both models (within a 
specified geometric tolerance) were identified through 
a comparison of the spatial geometry between the 
test and analysis graphics models. This resulted in 
streamlined analytical models which could then be 
compared in a direct manner with the test mod- 
els, thereby allowing a one-to-one comparison of grid 
point deformations. 

D-shape spar . Because the grid points along 
the leading edge of the D-shape spar were not ex- 
cited in the vertical direction, they appeared effec- 
tively constrained in this direction for each mode. 
Correspondingly, this lack of vertical deformation did 
not lend itself well to comparison with the associated 
analytical grid point deformations, which were free 
to move in the vertical direction, and thus degraded 
test-analysis correlation. To compensate for this, the 
vertical deformations of the associated analysis grid 
points were effectively negated for the cross-mapped 
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geometry model in the calculation of the MAC ma- 
trix for VD1. The analytical data sets containing the 
resultant deformations for the finite-element models 
were modified such that zero vertical deformation ex- 
isted for the selected grid points. This permitted the 
analysis model to mimic the lack of vertical defor- 
mation at the prescribed grid points while being free 
from the constraining effects associated with an ap- 
plied boundary condition. The resulting test-analysis 
MAC matrix for the D-shape spar is presented in fig- 
ure 41. The matrix is shown in a three-dimensional 
perspective view as well as a two-dimensional view to 
aid in the distinction of the constituent matrix mem- 
bers. Excellent agreement was obtained for the first 
vertical bending, first lateral bending, and first tor- 
sion modes, which correspond to the (1,1), (2,2), and 
(5,5) test-analysis matrix members. The associated 
MAC values for these pairs were in the range from 0.8 
to 1.0, as shown in table 21. The agreement between 
the test-analysis second vertical bending and second 
lateral bending modes was relatively fair, with MAC 
values of 0.565 and 0.484, respectively. The mag- 
nitude of these values seems somewhat low and is 
thought to be associated with slight differences in 
deformation characteristics for each mode. This is 
particularly evident in the comparison of the second 
lateral bending modes, where the analytical mode 
exhibits some torsional deformation which is not ap- 
parent in the test mode shape. Nevertheless, these 
modes have the characteristics of the second lateral 
bending shape and are judged to correspond to such. 
The additional modes, which correspond to the shell- 
type modes, occupy the remainder of the matrix and 
have MAC values from 0.0 to 0.8. It is noted that 
the MAC members corresponding to the test-analysis 
(7,1) and (1,6) matrix members have MAC values in 
the range from 0.6 to 0.8. In the former pair, the sev- 
enth test mode actually does resemble the first anal- 
ysis mode (first vertical bending). However, geomet- 
ric aliasing is promoting this elevated MAC value. In 
the latter pair, the first test mode (first vertical bend- 
ing) is similar to the sixth analytical mode, which is a 
shell-type mode. This similarity occurs as a result of 
the modification of the vertical deformation in ana- 
lytical grid points 7-12. Therefore, the cross-mapped 
model and the test model do appear similar in this 
mode, but this is clearly a consequence of compen- 
sating for the apparent immobilities of the test grid 
points and does not constitute a true agreement of 
mode shapes. 

The overall agreement for the D-shape tube fre- 
quencies is good, generally ranging from 1.9 percent 
difference (first vertical bending case) to 13.8 per- 
cent difference (second vertical bending case). The 


average difference between test and analysis 
frequencies is approximately 6.0 percent. 

Square spars . The interpretation of the MAC 
matrices for the square spars was straightforward 
in that all test and analysis grid points already 
possessed the same degrees of freedom; therefore, no 
additional adjustments needed to be performed. 

The test-analysis MAC matrix for VSl is shown 
in figure 42. Excellent agreement is exhibited for 
the first three global bending modes (first vertical, 
first lateral, and second vertical bending), with as- 
sociated MAC values ranging from 0.8 to 1.0, as 
listed in table 22. The second lateral bending mode, 
however, agrees poorly with its analytical counter- 
part, having a calculated MAC value of 0.243. A 
comparison of the test and cross-mapped analysis 
mode shapes, shown in figure 43, reveals that the 
two modes generally have the same shape. However, 
the cross-mapped analysis model is lacking the tor- 
sional deformation of the grid points farthest from 
the origin which is present in the test. This ap- 
pears to be the only significant difference between 
the two mode shapes, and it was therefore concluded 
that the test and analysis mode shapes do corre- 
spond, although some uncertainty exists. The first 
torsion mode from the test was found to agree fairly 
well with its analytical counterpart and had a MAC 
value of 0.573. Although the agreement of test- 
analysis data for VSl was within 9.5 percent differ- 
ence (second lateral bending), the average frequency 
difference was approximately 7.3 percent, with the 
best-case (first lateral bending) showing a 5.8 percent 
difference between test and analysis data. 

The test-analysis MAC matrix for VS2 is shown 
in figure 44. Excellent agreement between test and 
analysis mode shapes is again evident for the first 
three global bending modes (first vertical bending, 
first lateral bending, and second vertical bending), 
and the associated MAC values are listed in table 23. 
The second lateral bending mode shape correlated 
poorly with that obtained from the analysis, having 
a MAC value of 0.347. The analytical torsion mode 
(1753.94 Hz) showed roughly equal MAC values with 
the 1550.65-Hz test mode (0.564) and the 1702.77-Hz 
test mode (0.531). By inspecting the plotted mode 
shapes further (fig. 45), the test mode at 1550.65 Hz 
was judged to be a closer match for the analytical 
mode in terms of overall mode shape characteristics. 

Overall, the agreement between test and analysis 
frequencies was within 11.6 percent, with the best 
agreement corresponding to the first vertical bending 
mode (6.9 percent). The average difference between 
test and analysis data was approximately 8.9 percent. 
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The MAC matrix for the VS3 spar exhibited the 
best overall correlation of the square spars and is 
shown in figure 46. The MAC values correspond- 
ing to the global test and analysis pairs are listed in 
table 24. Excellent correlation was again obtained 
for the first three global bending modes, with the 
calculated MAC values ranging from 0.861 to 0.937. 
Comparison of the test and analysis shapes for sec- 
ond lateral bending mode (fig. 47) involved two test 
modes which appeared to share similar characteris- 
tics with the analytical mode shape. The test mode 
at 1423.28 Hz showed favorable correlation with the 
analytical mode at 1555.08 Hz in terms of calcu- 
lated MAC value (0.504), However, although the 
test mode corresponding to 1459.78 Hz shared simi- 
lar shape characteristics with its analytical counter- 
part, the associated MAC value for this analysis and 
test pair was only 0.011. This would seem to in- 
dicate a total lack of correlation for the two mode 
shapes, although some similarities are apparent from 
figure 47. Based on the relative MAC values, the 
mode at 1423.28 Hz better approximated the ana- 
lytical second lateral bending deformation. The re- 
spective test and analysis torsion modes correlated 
favorably, with a MAC value of 0.645. 

The agreement between the global test and anal- 
ysis frequencies was generally very good for the 
VS3 spar. The best agreement was found for the 
first lateral bending mode, showing a difference of 
only 2.3 percent. Although some uncertainty was 
observed regarding the test second lateral bending 
mode, the correlation between the test mode at 
1423.28 Hz and the analytical mode exhibited higher 
MAC values and agreed within 8.5 percent. The av- 
erage difference between test and analysis data was 
4.4 percent. 

Elliptical spars . Similar to the VD1 spar, the 
cross-mapped analysis models for the elliptical spars 
were modified to compensate for the lack of lateral 
deformation of the corresponding test grid points. In 
this case, the lateral deformations for the upper row 
of grid points (grid points 1-6 in fig. 9(c)) were zeroed 
out, as these grid points did not experience lateral 
excitations in the test. The resultant test- analysis 
MAC matrix for the VE1 spar is shown in figure 48. 
Excellent correlation was obtained for the five global 
modes. These modes had MAC values from 0.820 
to 0.980, as shown in table 25. The analytical mode 
shape corresponding to third vertical bending showed 
good correlation with the test mode at 1865.52 Hz, 
with a MAC value of 0.677. 

The data correlation for the VE1 spar exhib- 
ited very good agreement with experiment, with 


results ranging from 0.10 percent (first torsion) to 
4.9 percent (second lateral bending) for the five 
modes. The average difference between test and 
analysis frequencies was about 2.5 percent. 

The calculated MAC matrix for the VE2 spar is 
shown in figure 49. Excellent correlation was again 
obtained for the first four global modes. The as- 
sociated MAC values, shown in table 26, ranged 
from 0.843 to 0.978 for these modes. The test- 
analysis correlation for the first torsion mode was 
good (MAC value of 0.753), whereas the analyt- 
ical third vertical bending mode exhibited rela- 
tively poor correlation with the experimental mode 
at 1739.82 Hz (MAC value of 0.394). 

In general, the test-analysis results for VE2 ex- 
hibited very good agreement. The agreement ranged 
from 1.9 percent difference (first vertical bending) 
to 7.7 percent (first lateral bending) for the five 
modes. The average difference between test and 
analysis data was approximately 4.8 percent. 

The MAC matrix calculated for the third ellip- 
tical spar, VE3, is shown in figure 50. The matrix 
closely resembled those obtained for VE1 and VE2. 
Excellent agreement was obtained for the five global 
modes, with the calculated MAC values ranging 
from 0.800 to 0.980 as listed in table 27. The test- 
analysis correlation between the analytical third ver- 
tical bending mode and the experimental mode at 
1838.58 Hz was good, with an associated MAC value 
of 0.712. 

The agreement between test and analysis frequen- 
cies was generally very good and ranged from 1.6 per- 
cent (first vertical bending) to 9.6 percent (first lat- 
eral bending) for the first five modes. The average 
difference between test and analytical frequencies was 
approximately 6.0 percent. 

Concluding Remarks 

Dynamic finite-element analyses of thin-walled 
composite tubular spars exhibiting extension-twist 
coupling have been conducted and compared with 
experimental data. The cross-sectional design of the 
spars corresponded to three types: square, D-shape, 
and elliptical. Because these designs were noncir- 
cular in nature, they were therefore prone to cross- 
sectional warping. Finite-element analysis models 
were constructed for each spar to account for indi- 
vidual variations in wall thickness and, thus, weight, 
which were present in each manufactured spar. Re- 
sults from normal mode and frequency analyses of the 
finite-element models were compared with those ob- 
tained from free- free dynamic tests of the fabricated 
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spars. Five global or nonshell-type modes were iden- 
tified for the D-shape spar in the frequency range 
from 0 to 2000 Hz: first vertical bending, first lat- 
eral bending, second vertical bending, first torsion, 
and second lateral bending. The agreement between 
test and analysis frequencies corresponding to the 
mode shapes was good, ranging from 1.9 percent 
(first vertical bending) to 13.8 percent (second lat- 
eral bending). The average difference for all modes 
was approximately 6.0 percent. 

Five global modes were similarly extracted for 
each of three square spars (VS1, VS2, and VS3) in 
the frequency range from 0 to 2000 Hz. The identified 
mode shapes corresponded to first vertical bending, 
first lateral bending, second vertical bending, second 
lateral bending, and first torsion. The overall agree- 
ment between test and analysis results for the first 
square spar, VS1, was within 9.5 percent, with the 
best case agreement of 5.8 percent (first lateral bend- 
ing). The average difference in the test and analy- 
sis frequencies was 7.3 percent. The test and analy- 
sis data correlation for VS2 ranged from 6.9 percent 
(first vertical bending) to 11.6 percent (first torsion), 
with the average difference calculated at 8.9 percent. 
The test-analysis correlation for the final square spar, 
VS3, agreed within 8.5 percent overall. The average 
difference was 4.4 percent. 

Three elliptical spars were similarly tested and 
analyzed. Five global modes were obtained for each 
spar in the frequency range from 0 to 2000 Hz: first 
vertical bending, first lateral bending, second verti- 
cal bending, second lateral bending, and first tor- 
sion. The agreement between test and analysis fre- 
quencies for the VE1 spar ranged from 0.10 percent 
(first torsion) to 4.9 percent (second lateral bending) 
for the five global modes. The average difference 
between test and analysis frequencies was approxi- 
mately 2.5 percent. Similar trends were observed for 
VE2, with agreement of test and analysis frequen- 
cies ranging from 1.9 percent difference (first verti- 
cal bending) to 7.7 percent (first lateral bending) for 


the first five global modes. An average difference 
of 4.8 percent was observed between test and anal- 
ysis frequencies. The agreement between test and 
analysis frequencies for VE3 ranged from 1.6 percent 
(first vertical bending) to 9.6 percent (first ■ lateral 
bending) for the first five modes. The average dif- 
ference between test and analytical frequencies was 
approximately 6.0 percent. 

These results suggest that the structural dynamic 
characteristics of thin-walled composite structures 
employing extension-twist coupling can be deter- 
mined within practical engineering accuracy. The 
fabrication process employed in this investigation has 
demonstrated the potential for achieving a high level 
of repeatability in the manufacturing of composite 
structures. Despite inconsistencies in the cured wall 
thicknesses of the elliptical spar subset, this method 
generally yielded spars of excellent overall geometric 
integrity. It is important to note that a certain level 
of experience with this type of fabrication process is 
required, as there are several variables unique to the 
rubber expansion method which directly control spar 
quality. Repeatability in the manufacturing process 
is essential for the reliable and efficient reproduction 
of composite structural components that will be in- 
corporated in future rotorcraft systems. This inves- 
tigation of tubular spars provides a necessary data 
base to promote understanding structural dynamic 
characteristics of elastically coupled composite rotor 
blade spars. Future studies of increased scope and 
complexity, including the investigation of the effect 
of large values of pretwist in elastically coupled struc- 
tures, should further contribute to this data base 
and provide a basis for the development of practi- 
cal composite tilt rotor blade designs incorporating 
extension-twist coupling. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
October 16, 1992 
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Appendix 

Convergence Study Results 

A convergence study of limited scope was initi- 
ated to verify the adequacy of the finite-element mod- 
els. Specifically, the number of spanwise elements 
was varied from 16 to 33 to determine the effects 
upon the natural frequencies, particularly those cor- 
responding to the higher modes. Because this in- 
vestigation was primarily concerned with the global 
mode shapes (bending and torsion) and not the shell- 
type modes, the number of elements used circumfer- 
entially to model the cross section was not changed. 
The percentage of degrees of freedom retained in the 
analyses via static condensation was held at approx- 
imately 13 percent for each model per the recom- 
mendations outlined in references 19 and 20. The 


convergence study was performed with the finite- 
element model of the square spar. The results of 
this study are shown in table 28 for four test cases 
designated by “Model j” where j = 16, 20, 25, 
and 33 spanwise elements. The results indicated 
that the frequencies corresponding to the model us- 
ing 25 spanwise elements differed by a maximum 
of 0.65 percent when compared with the model using 
33 elements in the spanwise direction. Furthermore, 
the average difference between these two models was 
approximately 0.39 percent. If the characteristic 
shape of the global bending and torsion modes is sim- 
ilar among all spar types (i.e., number of elements 
necessary to properly define bending and torsional 
mode shapes), the conclusion can be made that the 
model using 25 elements in the spanwise direction 
was sufficient for this investigation. 
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Table 1. Material Properties of T650-42/ERLX 1925-2 
Graphite-Epoxy Cloth 


Eu, &22, lb/in 2 
G 12 , lb/in 2 . . 

^12 

p, lb/in 3 . . . 
^ply> • • • 


11.65 x 10 6 
8.20 x 10 5 
. . . 0.05 
. . 0.056 

. . 0.0075 


Table 2. Spar Wall Thicknesses and Weight Data 


Spar 

Average wall 
thickness, 
in. 

Maximum 

deviation, 

percent 

Measured 

weight, 

lb 

VDl 

0.0300 

0.0 

0.1896 

VD2 

0.0307 

2.3 

0.1962 

VD3 

0.0300 

0.0 

0.1929 

VS1 

0.0304 

1.3 

0.1929 

VS2 

0.0304 

1.3 

0.1940 

VS3 

0.0300 

0.0 

0.1918 

VE1 

0.0348 

13.8 

0.1631 

VE2 

0.0322 

6.8 

0.1598 

VE3 

0.0328 

8.5 

0.1731 


Table 3. Mounting Block Types and Weights 


Block type 

Weight, lb 

Accelerometer mount (D-shape, squares) 

0.0128 

Accelerometer mount (ellipses) 

0.0053 

Impact mount (leading edge of ellipses) 

0.0016 

Impact mount (upper surface of ellipses) 

0.0016 





Table 4. VD1 Test Node Locations 


Table 5. VS1 Test Node Locations 


Test 

node 

Global coordinates 

{x,y,z) 

1 

(0.63,-1.23,0.46) 

2 

(4.50,-1.23,0.46) 

3 

(10.13,-1.23,0.46) 

4 

(14.13,-1.23,0.46) 

5 

(19.50,-1.23,0.46) 

6 

(23.63,-1.23,0.46) 

7 

(0.63,0.38,0.00) 

8 

(4.50,0.38,0.00) 

9 

(10.13,0.38,0.00) 

10 

(14.13,0.38,0.00) 

11 

(19.50,0.38,0.00) 

12 

(23.63,0.38,0.00) 


Table 6. VS2 Test Node Locations 


Test 

node 

Global coordinates 

{x,y, z) 

1 

(0.63,-0.59,0.59) 

2 

(4.50,-0.59,0.59) 

3 

(9.25,-0.59,0.59) 

4 

(15.00,-0.59,0.59) 

5 

(19.50,-0.59,0.59) 

6 

(23.00,-0.59,0.59) 

7 

(0.63,0.59,0.59) 

8 

(4.50,0.59,0.59) 

9 

(9.25,0.59,0.59) 

10 

(15.00,0.59,0.59) 

11 

(19.50,0.59,0.59) 

12 

(23.00,0.59,0.59) 


Test 

node 

Global coordinates 

(x,y,z) 

i 

(0.63,-0.59,0.59) 

2 

(4.50,-0.59,0.59) 

3 

(9.38, -0.59, 0.59) 

4 

(14.13,-0.59,0.59) 

5 

(19.50,-0.59,0.59) 

6 

(22.88,-0.59,0.59) 

7 

(0.63,0.59,0.59) 

8 

(4.50, 0.59, 0.59) 

9 

(9.38,0.59,0.59) 

10 

(14.13,0.59,0.59). 

11 

(19.50,0.59,0.59) 

12 

(22.88, 0.59, 0.59) 


Table 7. VS3 Test Node Locations 


Test 

node 

Global coordinates 

{x,y,z) 

1 

(1.00,-0.59,0.59) 

2 

(4.50,-0.59,0.59) 

3 

(10.13,-0.59,0.59) 

4 

(13.75,-0.59,0.59) 

5 

(19.50,-0.59,0.59) 

6 

(23.00,-0.59,0.59) 

7 

(1.00,0.59,0.59) 

8 

(4.50,0.59,0.59) 

9 

(10.13,0.59,0.59) 

10 

(13.75,0.59,0.59) 

11 

(19.50,0.59,0.59) 

12 

(23.00,0.59,0.59) 
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Table 8. VE1 Test Node Locations 


Table 9. VE2 Test Node Locations 


Global coordinates 
(x,y, z) 


(1.00,-0.31,0.45) 

(4.50,-0.31,0.45) 

(9.50,-0.31,0.45) 

(14.13,-0.31,0.45) 

(19.50,-0.31,0.45) 

(23.00,-0.31,0.45) 

(1.00,-0.63,0.00) 

(4.50,-0.63,0.00) 

(9.50,-0.63,0.00) 

(14.13,-0.63,0.00) 

(19.50,-0.63,0.00) 

(23.00,-0.63,0.00) 


Test 

node 

Global coordinates 

(x,y,z) 

1 

(1.00,-0.31,0.45) 

2 

(4.75,-0.31,0.45) 

3 

(9.50,-0.31,0.45) 

4 

(14.50,-0.31,0.45) 

5 

(19.75, —0.31, 0.45) 

6 

(23.00,-0.31,0.45) 

7 

(1.00,-0.63,0.00) 

8 

(4.75,-0.63,0.00) 

9 

(9.50,-0.63,0.00) 

10 

(14.50,-0.63,0.00) 

11 

(19.75,-0.63,0.00) 

12 

(23.00, 1 -0.63, 0.00) 


Table 10. VE3 Test Node Locations 


Global coordinates 
(z, y, z) 

(1.00,-0.31,0.45) 

(4.50,-0.31,0.45) 

(9.50,-0.31,0.45) 

(14.50,-0.31,0.45) 

(19.50,-0.31,0.45) 

(23.00,-0.31,0.45) 

(1.00,-0.63,0.00) 

(4.50,-0.63,0.00) 

(9.50,-0.63,0.00) 

(14.50,-0.63,0.00) 

(19.50,-0.63,0.00) 

(23.00,-0.63,0.00) 











Table 11. Experimental Modal Parameters From Real Domain 
Poly reference Curve Fit for VD1 


Mode 

shape 

Frequency, 

Hz 

Damping, 

percent 

Amplitude 

Phase, 

deg 

MCF 

1 

“394.38 

0.372 

7.82 x 10 4 

-90.0 

0.999 

2 

“575.05 

0.356 

5.65 x 10 2 

90.0 

1.000 

3 

“904.19 

0.515 

1.69 x 10 4 

-90.0 

0.998 

4 

1051.48 

1.189 

1.56 x 10 5 

-90.0 

0.984 

5 

“1175.92 

0.444 

7.98 x 10 4 

-90.0 

0.999 

6 

1275.25 

0.448 

2.61 x 10 4 

-90.0 

0.997 

7 

1385.92 

0.445 

1.66 x 10 5 

-90.0 

1.000 

8 

1482.19 

0.383 

1.29 x 10 5 

-90.0 

0.994 

9 

1545.31 

0.381 

3.53 x 10 4 

-90.0 

0.979 

10 

“1576.74 

0.541 

8.51 x 10 4 

-90.0 

0.989 

11 

1769.89 

0.549 

4.23 x 10 4 

-90.0 

0.962 

12 

1897.91 

0.536 

2.70 x 10 3 

90.0 

0.979 

13 

1955.92 

0.533 

1.33 x 10 4 

90.0 

0.992 


“Global mode. 


Table 12. Experimental Modal Parameters From Real Domain 
Polyreference Curve Fit for VSl 


Mode 

shape 

Frequency, 

Hz 

Damping, 

percent 

Amplitude 

Phase, 

deg 

MCF 

1 

“464.05 

0.360 

7.27 x 10 4 

-90.0 

1.000 

2 

“486.27 

0.462 

2.05 x 10 4 

-90.0 

0.999 

3 

994.51 

0.665 

1.11 x 10 4 

90.0 

0.996 

4 

“1260.92 

0.437 

1.65 x 10 5 

-90.0 

0.999 

5 

1323.75 

0.287 

2.96 x 10 3 

-90.0 

0.985 

6 

1374.06 

0.230 

4.63 x 10 2 

-90.0 

0.998 

7 

“1412.16 

0.469 

1.11 x 10 5 

-90.0 

1.000 

8 

1435.64 

0.268 

2.32 x 10 5 

-90.0 

0.999 

9 

1591.98 

0.358 

8.80 x 10 3 

-90.0 

0.994 

10 

“1615.24 

0.394 

5.01 x 10 3 

-90.0 

0.999 


“Global mode. 
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Table 13. Experimental Modal Parameters Prom Real Domain 
Polyreference Curve Fit for VS2 


Mode 

shape 

Frequency, 

Hz 

Damping, 

percent 

Amplitude 

Phase, 

deg 

MCF 

1 

“465.12 

0.334 

8.82 

X 

10 5 

-90.0 

1.000 

2 

“478.51 

0.423 

2.01 

X 

10 5 

-90.0 

0.999 

3 

929.53 

0.346 

1.23 

X 

10 4 

-90.0 

1.000 

4 

“1240.09 

0.464 

2.57 

X 

10 6 

-90.0 

0.999 

5 

1267.61 

0.311 

5.01 

X 

10 3 

-90.0 

0.995 

6 

1330.98 

0.375 

2.27 

X 

10 6 

-90.0 

0.998 

7 

“1390.62 

0.415 

4.94 

X 

10 3 

90.0 

0.991 

8 

1478.98 

0.350 

6.14 

X 

10 2 

-90.0 

0.993 

9 

“1550.65 

0.459 

4.94 

X 

10 5 

-90.0 

0.997 

10 

1582.19 

0.295 

4.86 

X 

10 3 

-90.0 

0.983 

11 

1702.77 

0.513 

1.10 

X 

10 6 

-90.0 

0.999 


a Global mode. 


Table 14. Experimental Modal Parameters Prom Real Domain 
Polyreference Curve Fit for VS3 


Mode 

shape 

Frequency, 

Hz 

Damping, 

percent 

Amplitude 

Phase, 

deg 

MCF 

i 

“485.33 

0.383 

8.47 

X 

10 & 

-90.0 

0.999 

2 

“503.38 

0.398 

2.34 

X 

10 5 

-90.0 

0.999 

3 

1005.10 

0.362 

1.36 

X 

10 5 

-90.0 

0.998 

4 

“1297.59 

0.489 

2.73 

X 

10 6 

-90.0 

1.000 

5 

1367.69 

0.396 

1.53 

X 

10 6 

-90.0 

0.999 

6 

1379.64 

0.375 

1.15 

X 

10 5 

-90.0 

0.978 

7 

“1423.28 

0.311 

2.69 

X 

10 5 

-90.0 

0.988 

8 

1459.78 

0.416 

9.16 

X 

10 3 

-90.0 

0.983 

9 

1608.94 

0.309 

2.98 

X 

10 5 

-90.0 

0.981 

10 

“1653.92 

0.443 

1.51 

X 

10 5 

-90.0 

0.995 


a Global mode. 





Table 15. Experimental Modal Parameters From Real Domain 
Polyreference Curve Fit for VE1 


Mode 

shape 

Frequency, 

Hz 

Damping, 

percent 

Amplitude 

Phase, 

deg 

MCF 

1 

“382.55 

0.391 

1.26 x 10 b 

-90.0 

1.000 

2 

“481.01 

0.377 

1.42 x 10 4 

-90.0 

0.999 

3 

“1034.56 

1.417 

2.10 x 10 6 

-90.0 

0.999 

4 

“1279.04 

0.662 

1.95 x 10 5 

-90.0 

0.998 

5 

“1613.54 

0.602 

6.07 x 10 5 

-90.0 

0.998 

6 

1844.53 

2.226 

3.35 x 10 5 

-90.0 

0.990 

7 

1865.52 

0.729 

4.78 x 10 6 

-90.0 

0.987 


“Global mode. 


Table 16. Experimental Modal Parameters From Real Domain 
Polyreference Curve Fit for VE2 


Mode 

shape 

Frequency, 

Hz 

Damping, 

percent 

Amplitude 

Phase, 

deg 

MCF 

1 

“375.25 

0.555 

1.23 x 10 b 

-90.0 

0.998 

2 

“464.87 

0.395 

1.42 x 10 4 

-90.0 

0.999 

3 

“1013.66 

1.197 

2.11 x 10 6 

-90.0 

0.999 

4 

“1244.94 

0.624 

1.46 x 10 5 

-90.0 

0.996 

5 

“1571.33 

0.673 

1.43 x 10 6 

-90.0 

0.996 

6 

1739.82 

0.750 

2.31 x 10 6 

-90.0 

0.987 


“Global mode. 


Table 17. Experimental Modal Parameters From Real Domain 
Polyreference Curve Fit for VE3 


Mode 

shape 

Frequency, 

Hz 

Damping, 

percent 

Amplitude 

Phase, 

deg 

MCF 

1 

“378.46 

0.594 

1.21 x 10 6 

-90.0 

0.994 

2 

“457.75 

0.487 

1.37 x 10 4 

-90.0 

0.988 

3 

“1023.89 

1.071 

2.37 x 10 6 

-90.0 

0.981 

4 

“1233.83 

0.600 

7.59 x 10 4 

-90.0 

0.995 

5 

“1514.81 

0.617 

1.91 x 10 5 

-90.0 

0.996 

6 

1838.58 

0.655 

4.31 x 10 6 

-90.0 

0.985 


“Global mode. 






Table 18. Calculated Frequencies 
From Finite-Element Analysis 
for D-Shape Spar 


Mode 

shape 

VD1 

frequency, 

Hz 

i 

“402.03 

2 

“596.24 

3 

“1048.71 

4 

1200.44 

5 

“1254.97 

6 

1469.82 

7 

“1647.74 

8 

1695.49 

9 

1731.15 

10 

1914.47 

11 

2113.61 

12 

2423.44 

13 

2608.75 


a Global mode. 




Table 19. Calculated Frequencies From Finite-Element 
Analysis for Square Spars 


Mode 

Frequency, Hz, for — 

shape 

VSl 

VS2 

VS3 

1 

a 499.82 

“499.82 

“499.08 

2 

“516.02 

“516.02 

“515.54 

3 

1169.83 

1169.83 

1159.44 

4 

“1345.26 

“1345.26 

“1343.50 

5 

1445.49 

1445.49 

1430.12 

6 

1543.60 

1543.60 

1524.07 

7 

“1559.64 

“1559.64 

“1555.08 

8 

1581.47 

1581.47 

1562.83 

9 

“1753.94 

“1753.94 

“1737.28 

10 

11 

1844.37 

1844.37 

2518.15 

1827.15 


a Global mode. 


Table 20. Calculated Frequencies From Finite-Element 
Analysis for Elliptical Spars 


Mode 

Frequency, Hz, for — 

shape 

VEl 

VE2 

VE3 

i 

“383.30 

“382.71 

“384.46 

2 

“504.56 

“503.69 

“506.16 

3 

“1063.54 

“1061.99 

“1066.41 

4 

“1344.87 

“1342.42 

“1349.48 

5 

“1615.38 

“1611.57 

“1624.76 

6 

“2143.16 

“2136.97 

“2153.91 

7 

2304.54 

2279.08 

2368.29 


“Global mode. 





Table 21. Calculated Modal Assurance Criteria Values 
for Global Test and Analysis Modes for VD1 


Mode shape 

Matrix element 

Calculated 
MAC value 

Test 

FEA 

First vertical bending 

i 

1 

0.971 

First lateral bending 

2 

2 

0.941 

Second vertical bending 

3 

3 

0.565 

First torsion 

5 

5 

0.889 

Second lateral bending 

10 

7 

0.484 


Table 22. Calculated Modal Assurance Criteria Values 
for Global Test and Analysis Modes for VS1 


Mode shape 

Matrix element 

Calculated 
MAC value 

Test 

FEA 

First vertical bending 

1 

1 

0.921 

First lateral bending 

2 

2 

0.923 

Second vertical bending 

4 

4 

0.862 

Second lateral bending 

7 

7 

0.243 

First torsion 

10 

9 

0.573 


Table 23. Calculated Modal Assurance Criteria Values 
for Global Test and Analysis Modes for VS2 


Mode shape 

Matrix element 

Calculated 
MAC value 

Test 

FEA 

First vertical bending 

i 

1 

0.945 

First lateral bending 

2 

2 

0.944 

Second vertical bending 

4 

4 

0.801 

Second lateral bending 

7 

7 

0.347 

First torsion 

9 

9 

0.564 


Table 24. Calculated Modal Assurance Criteria Values 
for Global Test and Analysis Modes for VS3 


Mode shape 

Matrix element 

Calculated 
MAC value 

Test 

FEA 

First vertical bending- 

1 

1 

0.936 

First lateral bending 

2 

2 

0.937 

Second vertical bending 

4 

4 

0.861 

Second lateral bending 

7 

7 

0.504 

First torsion 

10 

9 

0.645 







Table 25. Calculated Modal Assurance Criteria Values 
for Global Test and Analysis Modes for VE1 


Mode shape 

Matrix element 

Calculated 
MAC value 

Test 

FEA 

First vertical bending 

1 

1 

0.980 

First lateral bending 

2 

2 

0.930 

Second vertical bending 

3 

3 

0.918 

Second lateral bending 

4 

4 

0.826 

First torsion 

5 

5 

0.820 

Third vertical bending 

7 

6 

0.677 


Table 26. Calculated Modal Assurance Criteria Values 
for Global Test and Analysis Modes for VE2 


Mode shape 

Matrix element 

Calculated 
MAC value 

Test 

FEA 

First vertical bending 

1 

1 

0.978 

First lateral bending 

2 

2 

0.922 

Second vertical bending 

3 

3 

0.926 

Second lateral bending 

4 

4 

0.843 

First torsion 

5 

5 

0.753 

Third vertical bending 

6 

6 

0.394 


Table 27. Calculated Modal Assurance Criteria Values 
for Global Test and Analysis Modes for VE3 


Mode shape 

Matrix element 

Calculated 
MAC value 

Test 

FEA 

First vertical bending 

i 

1 

0.980 

First lateral bending 

2 

2 

0.899 

Second vertical bending 

3 

3 

0.930 

Second lateral bending 

4 

4 

0.800 

First torsion 

5 

5 

0.822 

Third vertical bending 

6 

6 

0.712 






Table 28. FEA Convergence Study Results for Square Spar 


-Mode 

shape 

Frequency, Hz, for — 

Model 

16 

Model 

20 

Model 

25 

Model 

33 

1 

501.50 

500.43 

499.82 

499.02 

2 

516.68 

516.68 

516.02 

515.15 

3 

1180.22 

1172.67 

1169.83 

1165.72 

4 

1359.89 

1349.09 

1345.26 

1336.50 

5 

1453.30 

1447.39 

1445.49 

1441.86 

6 

1544.74 

1542.56 

1543.60 

1538.11 

7 

1573.29 

1562.94 

1559.64 

1550.91 

8 

1580.67 

1581.17 

1581.47 

1575.51 

9 

1772.39 

1759.15 

1753.94 

1746.49 

10 

1854.56 

1846.46 

1844.37 

1833.26 




Square 


D-shape 


Ellipse 


Figure 1. Cross-sectional designs for composite spars. Dimensions are in inches. 
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Figure 2. Initial composite mandrel wrap and silicon rubber tool. 
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Figure 3. Initial wrapping of composite material. 
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Figure 4. Continuation of wrap procedure using subsequent composite layers. 
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Figure 5. Female aluminum tool cavity and cover plate. 



Figure 6. Cured tubular spars. 
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Figure 7. Vibration test setup. 



Figure 8. Triaxial accelerometer configuration. 
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(a) D-shape spar. 


(b) Square spar. 



Figure 9. Degrees of freedom for test spars. 



Figure 10. Data acquisition system. 
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grid points 
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dynamic degrees of freedom 


Figure 11. Finite-element model of D-shape spar. 



Figure 12. Finite-element model of square spar. 
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Figure 13. Finite-element model of elliptical spar. 
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Figure 16. 1 Z_ driving point measurement for VD1. 







(a) First vertical bending. 


(b) First lateral bending. 










(c) Second vertical bending. 


(d) First torsion. 




(e) Second lateral bending. 

Figure 17. Global mode shapes for VD1. 
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Figure 18. Multivariate mode indicator function for VS1. 
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Figure 19. Multivariate mode indicator function for VS2. 
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Figure 21. Frequency response function for VS1. 
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Figure 22. 1 Z- driving point measurement for VS1. 







(a) First vertical bending. (b) First lateral bending. 



(c) Second vertical bending. (d) Second lateral bending. 






iZ 



(e) First torsion. 

Figure 23. Global mode shapes for VSl. 
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Figure 25. 1 Z- driving point measurement for VS2. 
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Figure 26. Frequency response function for VS3. 
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Figure 27. 1 Z- driving point measurement for VS3. 


41 








MMIF 


MMIF 





MMIF 



Phase, 

deg 




0 500 1000 1500 2000 

Frequency, Hz 

Figure 31. Frequency response function for VE1. 







Frequency, Hz 

Figure 32. 1 Z- driving point measurement for VE1. 


44 







(c) Second vertical bending. (d) Second lateral bending. 





(e) First torsion. (f) Third vertical bending. 

Figure 33. Global mode shapes for VE1. 
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Figure 35. 1 Z- driving point measurement for VE2. 
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Figure 36. Frequency response function for VE3. 
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Figure 37. 1 Z- driving point measurement for VE3. 









(a) First vertical bending. (b) First lateral bending. 



(c) Second vertical bending. (d) First torsion. 



(e) Second lateral bending. 


Figure 38. Global analytical mode shapes for VD1. 
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(b) First lateral bending. 



% 


x J, 



% 


'J, 


(c) Second vertical bending. 


(d) Second lateral bending. 
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(e) First torsion. 

Figure 39. Global analytical mode shapes for VS1. 
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(a) First vertical bending. 


(b) First lateral bending. 



(c) Second vertical bending. 



(d) Second lateral bending. 



(e) First torsion. (f) Third vertical bending. 

Figure 40. Global analytical mode shapes for VE1. 
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Figure 41 . Modal assurance criteria matrix for VD1 . 



Figure 42. Modal assurance criteria matrix for VS1. 
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Frequency = 1559.64 Hz 



(a) Cross-mapped VS 1 finite-element model. 


(b) VS1 test model. 


Figure 43. Test-analysis second lateral bending mode shapes for VS1. 
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Figure 44. Modal assurance criteria matrix for VS2. 


ORIGINAL PAGE 
COLOR PHOTOGRAPH 


52 





J 

x-«— Jy 

(a) Cross-mapped VS2 finite-element model. 
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(b) VS2 test models. 
Figure 45. Test-analysis torsion mode shapes for VS2. 



Figure 46. Modal assurance criteria matrix for VS3. 
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Frequency = 1555.08 Hz 

(a) Cross-mapped VS3 finite -element model. 

Figure 47. Test-analysis second 



(b) VS3 test models, 
lateral bending mode shapes for VS3. 



Figure 48. Modal assurance criteria matrix for VE1. 
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Figure 49. Modal assurance criteria matrix for VE2. 
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Figure 50. Modal assurance criteria matrix for VE3. 
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